N90 



-23024 


Modern CACSD using the Robust-Control 

Toolbox 


Richard Y. Chiang and Michael G. Safonov 
Department of Electrical Engineering - Systems 
University of Southern California 
Los Angeles, CA. 90089-0781 

i 


Abstract 

The Robust-Control Toolbox [1] is a collection of 40 “M-files” which ex- 
tend the capability of PC/PRO-MATLAB to do modem multivariable robust 
control system design. Included are robust analysis tools like singular values 
and structured singular values, robust synthesis tools like continuous/discrete 
H 2 /H°° synthesis and LQG Loop Transfer Recovery methods and a variety of 
robust model reduction tools such as Hankel approximation, balanced trunca- 
tion and balanced stochastic truncation, etc. 

In this paper, we will describe the capabilities of our toolbox and illustrate 
them with examples to show how easily they can be used in practice. Examples 
include structured singular value analysis, H°° loop-shaping and large space 
structure model reduction. 


1 Introduction 

The fundamental issue in robust control theory - to find a stabilizing controller that 
achieves feedback performance despite the plant uncertainty, is still the same issue 
addressed by the classical 1930’s feedback theory of Black, Bode and Nyquist (ref. 
Fig. 1.1). Modern robust control theory has resolved many of the issues concerning 
the “gap” between the theory and practice that had grown to troublesome proportions 
in the 1970’s. One key to bridging the “gap” has been the singular value Bode plat. 
Recent progress in Structured Singular Value (SSV), H°° optimal control theory and 
the model reduction techniques utilizing singular values have made the modern robust 
control theory highly practical. 

The inavailability of quality software implementing the techniques of robust con- 
trol theory has, until very recently, significantly limited the access of both researchers 
and engineering practitioners to these techniques. 
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Fig. 1.1 Robust Control Problem. 

Our Robust-Control Toolbox implements the most up-to-date robust control theory 
like Perron SSV, optimal descriptor 2-Riccati H°° formulae, LQG/LTR and singular 
value based model reduction techniques, . . . , etc. The toolbox consists of a library 
of 40 functions which extend the capabilities of PC I PRO - MAT LAB™ and the 
PC / PRO M AT LAB Control Toolbox. The toolbox itself represents four man years 
research work done at USC. 

These 40-functions can be catalogued into 3 major areas: 

• Robust Analysis 

— Singular Values 

— Characteristic Gain Lod 

- Structured Singular Values 

• Robust Synthesis 

— LQG/LTR, Frequency- Weighted LQG 

- H 2 , H°° 

• Robust Model Reduction 

- Optimal Descriptor Hankel (with Additive Error Bound) 

- Schur Balanced Truncation (with Additive Error Bound) 
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- Schur Balanced Stochastic Truncation (with Multiplicative Error Bound) 

In this paper, we will highlight the most important functions in the toolbox, 
demonstrate how easily they can be used and show what kind of results can be 
achieved with practical problems. As for the details of the modern robust control 
theory, they can be found in [7], [2] and the references therein. 

2 Robust Analysis 

The objective of robust analysis is to find a proper measure of the multivariable 
stability margin (MSM) against uncertainty. Uncertainty may take many forms, 
but among the most significant ones are noise/disturbance signals, transfer function 
modeling errors and unmodeled nonlinear dynamics, etc. Uncertainty in any form is 
no doubt the major issue in most control system designs. 

Several tools to measure MSM are available: [1], [8] 

• Singular values (Safonov, 1977; Doyle, 1978) 

• Perron eigenvalues (Safonov, 1982) 

• Diagonal scaling via nonlinear programming (Doyle, 1982; Tekawy et al., 1989) 
Let’s define the MSM first: 

Definition 1 Multivariable Stability Margin (MSM)K m , fi(-)" 1 

K m (Tyu) = KTyu)' 1 = inf{<f(A)| det(I - T^A) = 0}. 

In other words, it’s the smallest <f(A) that can make the determinant (I — TyuA) 
singular (or the closed-loop system unstable). See Fig. 1.1. 

A theorem summarizes the whole MSM idea: 

Theorem 1 The system is stable for all stable A* with || A^l^ < 1, if the MSM 

K m cr m )> i. 

Unfortunately, exact computation of K m (or /i -1 ) would require solution of a non- 
con vex optimization problem and is therefore impractical. Fortunately, computable 
upper bounds on K m are available, viz., 

K- l (T^) = KTyu) < mf ||DT vu ir 1 || 00 < mf ||0aMT yu )ir 1 || oo < | I^WIU 

where V := {diag(dil , . . . ,d n I)\di > 0}. 

Then using these upper bounds (some may be more conservative than others), one 
can assure that the system is stable against the norm-bounded uncertainty || A < 
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A comparison of the available upper bounds reveals that some are much easier to 
compute than others. See table 2.1. 

Table 2.1 


Method 

Property 

Computation 

Reference 

Optimal Diagonal 
Scaling 

n = 3, exact K m 
n > 3, 3 15 % gap 

demanding 

Doyle, 1982 
Tekawy et al., 1989 

Perron Eigenvector 
Diagonal Scaling 

very close to optimal 
diagonal scaling 

easy 

Safonov, 1982 

Singular Value 

can be very 
conservative 

easy 

Safonov, 1977 
Doyle, 1978 


Let’s see the following example. 

Example: Given a system G(s) with multiplicative uncertianty at its input. Find 
the MSM. 

JL 

. * 4-4 “ I " *+8 *+8 

Theorem 1 implies ||A|U < (||G(J + G)" 1 || 00 )“ 1 . 



KaiVSec 


Fig. 2.2 Singular Value vs. K m (upper bound). 

This example reveals that the singular value upper bound is too conservative in 
“robust analysis”. Whereas, Perron SSV is much simpler to compute than diagonally 
scaled nonlinear programming fi. 


3 Robust Synthesis 

Classical control system designers often do “loop-shaping” to meet design specifi- 
cations. So do modern robust control system designers. “Loop-s haping ” for mul- 



















tivariable systems is done via the singular- value Bode plot. However, to shape the 
loop transfer function L(s) is nothing but to shape the sensitivity function S(s) = 
(I + L(s ))~ x and the complementary sensitivity function T(s) = L(s)(I + L(s)) 1 . 
See Fig. 3.1 



There are several loop-shaping methods available in the Robust- Control Toolbox 
(see Table 3.1), but H°° is one of our favorites. 



Disadvantages 


o need full-state feedback 
o need accurate model 
o possibly many iterations 


o no stability margin guaranteed 
o need accurate model 
o possibly many iterations 


o high gain controller 
o possibly many iterations 
o design focus on one point 


o possibly many iterations 


Example: Classical loop shaping vs. H°° for 2nd order low-damped system. 

Given a plant Gr($) which is 2nd order with damping 0.05 at 20 rad/sec, find a 
controller to meet frequency response Bode plot (see Fig. 3.2) 
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PLANT OPEN LOOP & SPECIFICATION 



Fig. 3.2 2nd order plant open loop (£ : 0.05,0.5) and the L(s ) spec. 

A classical design might be decomposed into the following: (see Fig. 3.3) 

Step 1: Rate feedback to improve damping. 

Step 2: Design high frequency (phase margin, BW, roll-off..). 

Step 3: Design low frequency (DC gain, disturbance rejection..). 

The classical result is shown in Fig. 3.4. Now, let’s see how H°° approaches the 
problem. 


H°° Problem Formulation 

We are solving the so-called H°° Small-Gain Problem ([3]) using the numerically 
robust “optimal” descriptor 2-Riccati formulae of Safonov, Limebeer and Chiang [4] 
[5]. 

H°° Small-Gain Problem: 

Given a plant P(s) (ref. Fig. S.5), find a stabilizing controller F(s) such that 
the closed-loop transfer function Ty^ Ul is internally stable and its infinity-norm is less 
than or equal to one. 

But what makes H°° work is its unique and remarkable “all-pass” property: 

At H°° optimal, the frequency response o/T WUl is all-pass and equal to one (i.e., 

Iir„„|| ^ i)i 

This means that designers can achieve EXACT frequency domain loop-shaping via 
suitable weighting strategies. For example, one may augment the plant with frequency 
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lead-lag 


plant 



Fig. 3.3 Classical loop-shaping block diagram. 


CLASSICAL LOOP-SHAPING 



Fig. 3.4 Classical loop- shaping. 
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Fig. 3.6 Weighti 
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dependent weights W\, W 2 and W 3 as shown in Fig. 3.6. Then, the Robust-Control 
Toolbox functions augtf.m and augss.m will perform the augmentation and create a 
state-space for function hinf.m to find an H°° controller. Of course, these frequency 
weighting functions have to be chosen so that a stabilizing solution satisfying the H°° 
norm constraint exists. 

In a typical application, either W^s) or W^a) would be absent, leading to weighted 
H°° costs of the forms 


min 

F(,) 


W 2 S ' 
W Z T 


< 1 

oo 


or 


min 

n*) 


WiS 

w 2 fs 


< 1 . 

oo 


In our example, the frequency domain spec, can be split into W\ and W 3 : 



(0.2s + l) 2 
^100(0.005s + l) 2 ’ 



40000 

s 2 


as shown in Fig. 3.7. 


SPECIFICATION — > WEIGHTING STRATEGY 



Fig. 3.7 Weighting strategy for 2nd order problem spec. 

The results are shown in Fig. 3.8 for different p’s. Clearly, in the limit (as p 
goes to 3.16) the cost function becomes “all-pass”. The parameter p of W\ is the 
only parameter on which we iterate for design; the Robust- Control Toolbox script-file 
hinfgama.m automates this iteration. 
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Fig. 3.8 H°° results for 2nd order system 
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3.1 H°° Software Execution and Sample Run 

To do an H°° control design with the Robust-Control Toolbox is relatively simple. 
Table 3.2.1 shows the complete user inputs for our sample problem. 

Table 3.2.1 

>> nug = [0 0 400]; dng = [1 2 400]; 

» [ag,bg,cg,dg\ = tf2ss(nug,dng)-, 

» sysg = [ag bg\cg dg ]; xg = 2; 

>> ti>l = [2.5e — 5 l.e — 2 1; i 

0.01 * [4.e - 2 4.e — 1 1]]; 

» w2 = [ ]; 

» u>3 = [1 0 0; 0 0 40000]; 

» [A,B1,B2,C1,C2,D11,D12,D21,D22] = augtf(sysg,xg,wl,w2,w3); 

>> hinf — 

Table 3.2.2 shows the output which appears on the screen for a successful run of 
hinf.m. 


Table 3.2.2 



« H-inf Optimal Control Synthesis >> 



Computing the 4-block H-inf optimal controller 
using the S-L-C loop-shifting/descriptor formulae 

Solving ior the H-inf controller F(s) using U(s) 

* 0 (default) - - - 

Solving riccati aquations and performing H-infinity 

existence tests: 


i. 

Is Oil small enough? 

OK 

2. 

Solving state-feedback (P) riccati . . . 



a. No Hamiltonian jv-axis roots? 

OK 


b. A-B2*F stable (P >■ 0)? 

OK 

3. 

Solving output -inject ion (S) riccati . - . 



a. No Hamiltonian jir-axis roots? 

OK 


b. A-G*C2 stable (S >= 0)? 

OK 

4. 

max eig(P*S) < 1 ? 

OK 

all tests passed — computing H-inf controller . 

. . 


DONE! ! ! 


\ , 
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4 Robust Model Reduction via Basis-Free Tech- 
niques 

In the design of controllers for complicated systems, model reduction arises in several 
places: 1). Plant model reduction, 2). Controller model reduction, 3). Simulation of 
large size problem. 

However, naive implementations of model reduction methods such as Rosenbrock’s 
stair-case algorithm, Moore’s balanced truncation, optimal Hankel approximation and 
balanced stochastic truncation, etc., can fail on even relatively simple problems due 
to numerical instability. 

The Robust-Control Toolbox implements the “basis-free” version of the latter three 
of the model reduction techniques, which are not only numerically robust but also 
tend to achieve the ultimate result for practical problems. In particular, they all 
possess the following special features: 

1. They bypass the ill-conditioned balancing transformation, so that they ran eas- 
ily deal with the “non-minimal” systems. 

2. They employ Schur decomposition to robustly compute the orthogonal bases 
for eigenspaces required in intermediate steps. 

These methods all enjoy attractive L°° error bounds - either an additive error 
bound or a multiplicative error bound. 

• Additive Methods: 

— Optimal descriptor Hankel MDA (ohklmr.m). 

— Schur balanced truncation (balmr.m, schmr.m). 

• Multiplicative Method: 

— Schur balanced stochastic truncation (reschmr.m). 


4.1 Robust Model Reduction Theorems 

For robust control system design, it is desirable that the reduced order model satisfies 
the conditions of one of the following two theorems (ref. [1] [2]). Otherwise, the 
controller design based on the “blind”reduced order plant model can be unstable ! 

Theorem 2 Additive Robustness Theorem; (see Fig. \.l) 

If < S.{G) for u> < u> r (with A a and G open loop stable), then the closed U 

loop system will be stable provided that the control bandwidth is less than u> r , where 
u> r := max{u> \ o(G(jw)) > o(A A (ju>))}. 
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TRUE PLANT G(s) 



Fig. 4.1 Additive modeling error. 

Theorem 3 Multiplicative Robustness JTheorem: (see Fig. 4-2) 

If < 1 for u> < u> r (with Am and G open loop stable), then the closed-loop 

system will be stable provided that the control bandwidth is less than u> r » where u> T := 
max{uj | 5"( Ajvf(yw)) < 1}. 


TRUE PLANT G(s) 



Fig. 4.2 Multiplicative modeling error. 


4.2 Examples of Model Reduction 

Example 1: Find a 3-state reduced order model for the transfer function 

0.05(s 7 + 801s 6 + 1024s 6 -I- 599s 4 + 451s 3 + 119s 2 + 49s + 5.55) 
~ s 7 + 12.6s 6 + 53.48s 6 + 90.94s 4 + 71.83s 3 + 27.22s 2 + 4.75s + 0.3 

with Hankel singular values of the phase matrix 


01 

02 

03 

04 

05 

ere 

cr 7 

0.9959 

0.9972 

0.9734 

0.7166 

0.5635 

0.0021 

0 
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Fig. 4.3 Schur BST-REM vs. Schur BT and Descriptor Hankel. 
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The results produced by the Robust-Control Toolbox functions ohklmr.m, schmr.m, 
reschmr.m, are shown in Fig. 4.3. Clearly, the Schur BST-REM method that keeps 
the reduced model staying inside a prescribed relative error bound produces the ul- 
timate result in model reduction. Note that <77 = 0 indicates only the “basis-free 
methods can handel the problem without numerical difficulty. 

Example 2: Model reduction for a large space structure [6] (see Fig. 4.4). 

Our design reequirement is to find a controller to track LOS loops in 300 Hz BW 
and to reduce plant disturbance response by a factor of 100. 

The Hankel singular values after the inner loops are closed indicate that the system 
is non-minimal. Therfore, only the “basis-free” methods such as — Schur BT (schmr.m) 
or Schur BST-REM (reschmr.m) can be used. Again, only the Schur BST-REM can 
match the original model up to a “robust frequency” which is high enough so that 
the required BW of 300 Hz can be satisfied (see Fig. 4.5 & Fig. 4.6). 
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Fig. 4.4 Large space structure. 
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LSS 4-tute Schur-Bt model vs. 116-scaae original 


LSS 4-state Schur-BT model vs. local error 
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Fig. 4.5 Model reduction using Schur Balanced Truncation. 


LSS 4-stale Schtv BST-REM model vs. 116-soae orig. 
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Fig. 4.6 Model reduction using Schur Balanced Stochastic Truncation. 
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5 Importance to the Control Community 


If the ultimate goal of research is to apply the theory to reality, then the contribution 
of Robust- Control Toolbox is clear: 

It provides a vital bridge between modern robust control theory and real control appli- 
cations. 

The following diagram (Fig. 5.1) shows how different groups of people with dif- 
ferent backgrounds can utilize the Robust- Control Toolbox to achieve their personal 
goals. For example classical control designers can use the toolbox to understand the 
theory or to apply on a design. Non-robust control researchers can study the code pro- 
vided by the toolbox together with the papers referenced therein to become familiar 
with the robust control theory. Students can use the toolbox either for robust control 
research or for realistic design studies. It seems that the Robust- Control Toolbox can 
serve people from a variety of backgrounds. 


Robust 

Control 

Toolbox 



Robust 




Fig. 5.1 The importance of the Robust- Control Toolbox. 
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6 Summary 

The Robust- Control Toolbox 

• Provides a bridge between modern robust control theory and the real-world 
applications. 

• Has the most up-to-date Robust-Control theories and algorithms. 

• Is in readable M-files, so it’s educational. 

• Contains the tools one needs to do robust control system design, analysis and 
model reduction. 

• Is direct, powerful and easy-to-use. 
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